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O ' 

ON ■ ABSTRACT 

^"'i We have evaluated a systematic effect on counts-in-cells analysis of deep, wide- field 

galaxy catalogues induced by the evolution of clustering within the survey volume. A 
multiplicative correction factor is explicitly presented, which can be applied after the 
higher order correlation functions have been extracted in the usual way, without taking 
into account the evolution. The general theory of this effect combined with the ansatz 
describing the non-linear evolution of clustering in simulations enables us to estimate 
the magnitude of the correction factor in different cosmologies. In a series of numerical 
calculations assuming an array of cold dark matter models, it is found that, as long as 
galaxies are unbiased tracers of underlying density field, the effect is relatively small 
(~ 10%) for the shallow surveys (z < 0.2), while it becomes significant (order of unity) 
in deep surveys (z ~ 1). Depending on the scales of interest, the required correction 
is comparable to or smaller than the expected errors of on-going wide-field galaxy 
surveys such as the SDSS and 2dF. Therefore at present, the effect has to be taken 
into account for high precision measurements at very small scales only, while in future 
deep surveys it amounts to a significant correction. 



Subject headings: cosmology: theory - dark matter - large-scale structure of universe - 
galaxies: distances and redshifts 



The Astrophysical Journal (Letters), in press 



- 2 - 



1. Introduction 

Cosmological observations are necessarily carried out on a null hypersurface or a light-cone. 
At low redshifts (z < 0.1), this can be regarded as to provide information on the constant-time 
hypersurface (z = 0) which is a quite conventional implicit approximation underlying cosmological 
studies using the galaxy redshift surveys. When the depth of the survey volume exceeds z ~ 0.1, 
however, this approximation breaks down, and one should simultaneously take account of the 
intrinsic evolution of galaxy clustering and the light-cone effect in addition to any other selection 
effect in interpreting the data. This is indeed the case for the on-going wide-field surveys of 
galaxies including 2dF (2-degree Field Survey) and SDSS (Sloan Digital Sky Survey). 

To our knowledge, the first quantitative consideration of the light-cone effect is made 
by Nakamura, Matsubara, k Suto (1998) who derived the systematic bias in the estimate of 
(3 ~ r2g 6 /6 from magnitude-limited surveys of galaxies combining the cosmological redshift 
distortion effect (Matsubara k Suto 1996) and the evolution of galaxy clustering within the 
survey volume. In this paper, we examine the light-cone effect on higher-order statistics of galaxy 
clustering, considering counts-in-cells analysis specifically. 

Let us consider first the higher-order statistics on the idealistic constant-time hypersurface. 
Denote the volume averaged N-th order correlation functions at a redshift z by £ N (R;z), 
where R is the comoving smoothing length, and introduce the normalized higher-order moments 
Sn(R;z) = £,n(R; z)/[^ 2 (R'i z )] N_1 - The hierarchical clustering ansatz states that Sn(R;z) is 
constant and independent of the scale R. This is a good approximation in nonlinear regimes, 
although small but definite scale-dependence is clearly detected from N-body experiments (Lahav 
et al. 1993; Suto 1993; Matsubara k Suto 1994; Suto k Matsubara 1994; Jing k Borner 1997). In 

r_ -I N-l 

addition, perturbation theory predicts that ^ N (R;z) evolves in proportion to ^ 2 {^ z ) > an d 
therefore Sn(R;z) is independent of time, i.e. it is constant with respect to z. 

The next section describes the general theory of the light cone effect on Sn(R;z) defined 
above. Using the ansatz by Jain, Mo, k White (1995; hereafter JMW), §3 evaluates the 
appropriate correction in an array of cold dark matter (CDM) models. Finally, §4 summarizes the 
results and discusses the implications for redshift surveys. 



2. Observing the higher-order moments on the light-cone 

It is difficult to estimate £jy (-^i z ) observationally since z is changing over the volume of 
galaxy sample. While in principle one could measure the iV-point functions on z = const surfaces, 
in practice this would result in a diminished volume, thus a significant increase of the errors. 
Instead it is more practical to extract the following N-th order correlation functions averaged over 
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the volumes on the light-cone: 



z 2 dz w(z)£ N (R; z) 

lax 

z dz w(z) 



(1) 



In the above expression, we assume that the observation is performed with the fixed solid angle, 
and the sampling cells for the analysis are placed randomly in z-coordinate with w(z) being their 
weighting function. If the cells were located randomly in the comoving coordinates, the volume 
element z 2 dz should have been replaced by dji(z; Oo, Xq) 2 c\dt / dz\dz (dA is the angular diameter 
distance; see, Nakamura et al. 1998) and thus the procedure itself becomes dependent on adopted 
values of S7o and Ao- In principle w(z) is an arbitrary function, and should be determined so as to 
maximize the signal-to-noise ratio given the selection function of individual observation. By z max 
we denote the redshift corresponding to the depth of the survey. For a volume-limited sample, for 
instance, it is natural to set w (z) = 1 and z max as the maximum redshift of the sample. Similarly 
we define the (observable) higher-order moments averaged over the light-cone as 



S N (R; < z n 



£ N (R;< z„ 



fiziR'i < ^max) 
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It is useful to introduce the function G(z) which describes the evolution of the averaged 
two-point correlation function: 

l(R-z) = G(zMR;0). (3) 

In linear regime, G(z) is equivalent to [D(z)/ D(0)] 2 where D(z) = D(z; Qq, Xq) is the linear growth 
rate. Although the above relation (^) does not exactly hold in the nonlinear regime, several 
approximation formulae are derived in the literature, which empirically describe the evolution by 
allowing G(z) depend on the scale R (see §3 for details). 

Once we accept the evolution law (||) , equation @ is explicitly written as 



S N (R; < z n 
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If ^max "C 1, the above expression is expanded in terms of z max as follows: 



Sn{R] < ^max 

3 



S N (R;0) + ^S' N (0)z n 
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(AT - 1)(JV - 2)S N (0)G'(0) 2 + — {N - 1)S' N (0)G'(0) + —S'^(0) 
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(4) 
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where S' N (0) denotes dS]y(R; z)/dz\ z= o and so on. The above expansion up to 0(z max ) is valid as 
long as the weighting function is well-approximated up to the same order: 

™(Zmax) = W(0) + w'(0)z max + -w" (0) Z max • ( 6 ) 

In other words, z max should be set to be smaller than the effective window size of u>(z max ). 

It is interesting to note that up to 0(z max ) equation (||) is independent of w(z), and that 

0(z 

max) term is determined only by Sjy(0) independently of G(z). Since S'j^(0) is expected to 
vanish in linear theory (Fry 1984; Goroff et al. 1986; Bernardeau 1992), and shown to be relatively 
small even in quasi- and fully- nonlinear regimes (Bouchet et al. 1992; Lahav et al. 1993; Colombi, 
Bouchet, & Hernquist 1995; Szapudi et al. 1997), equation (||) implies that the light-cone effect is 
very small for 2dF and SDSS galaxy redshift surveys (z max < 0.2). It should be noted, however, 
that if galaxies are biased relative to the mass density field, S' N (0) may not be necessarily small. 
So any signal proportional to -z max provides a clear indication of the time-dependent biasing of 
galaxies (see e.g., Fry 1996; Mo & White 1996; Mo, Jing k White 1997). 



3. Evaluating the light-cone effect: an example 

In order to evaluate the effect of observational average on the light-cone, we assume that 
Sn(R] z) does not evolve with z, i.e., Sn(R; z) = Sn(R; 0). As described above, this is a reasonable 
approximation as long as galaxies are unbiased tracers of underlying density field. If we introduce 
the measure of the light-cone effect: 

A N (R, < z max ) = SN{m -1, (7) 

equations (H) and (0) with Sn(R;z) = Sn(R;0) reduce to 

Ajv(w) = -^(N -1)(N- 2)G'(0) 2 4 ax + 0(z max ). (8) 

Note that (1 + Ajv) can be regarded as a correction factor as well, if one measures the Sn's without 
considering the evolution of clustering. This constitutes a simple and practical method, which we 
propose for future measurements, when compensation for the light cone effect is required. 

To evaluate equation (||), we need a model for G(z). For this purpose, we adopt the ansatz 
originally put forward by Hamilton et al. (1991) and improved later by Peacock & Dodds (1994, 
1996) and JMW. To be specific, we apply the fitting formula by JMW which relates the evolved 
two-point correlation function £#(i2; z) with its linear counterpart £l(-Ro! z) as follows: 

l E {R;z) = B(n)F[l L (Ro;z)/B(n)], (9) 
. . x + 0.45x 2 - 0.02x 5 + 0.05x 6 . , 

F{X) = 5 ■ (10) 

v ' 1 + 0.02x 3 + 0.003x 9 / 2 v ; 
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In the above equations, n denotes a power-law index of the power spectrum, Rq = [l+£ E (z, i?)] 1//3 i2, 
and B(n) = [(3 + re)/3] ' 8 . JMW show that the above formula works reasonably well even for 
CDM models by replacing n by the effective spectral index evaluated at the scale which is just 
entering nonlinear regime. In general the resulting n depends on z, which we neglect below for 
simplicity; for galaxy surveys which we are primarily interested in, the z-dependence of n near 
z = is expected to be very small. 

The inverse of equation (||) is formally written as ^ l (Rq;z) = B(n)F~ 1 [^ E (R; z) / B{n)\ and 
JMW's empirical fit to F~ 1 (y) is 



F-\y) 



1 + 0.036y iyd + 0.000V 
1 + 1.75?/ - 0.0015y 3 - 63 + 0.028y 4 



1/3 



(11) 



Then £e(R; z) is expressed explicitly in terms of (, E (R; 0): 



H E (R; z) = B(n)F 



D 2 (0) 



F-^ E (R,0)/B(n)} 



(12) 



Let us introduce a parameter a(R) = F^ 1 £ E (R,0)/B(n) which characterizes the variance 
on a scale R at z = 0, and thus depends on f^o and Ao through the shape of the fluctuation 
spectrum. Then the scale-dependent evolution factor G(z) = G(R; z) in equation (0) is given by 



G(R;z) 



1 



£ E {R;0) ' F(aY 
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For the convenience of z-expansion, we calculate the derivatives of the above quantity at z = 0: 



dG(R;z) 



dz 

d 2 G(R; z) 
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(14) 
(15) 



where /o = din D/da\ z= o, qo = £Iq/2 — \$. The above expressions indicate how the light-cone 
effect depends on Q,q and Ao at z max -C 1. Note, that they are involved in 0(z 2 a£ix ) term and thus 
do not contribute significantly at small z. 



4. Results and conclusions 

Using equations (f|) and ( |i~3| ) and assuming Sn(z) = Sn{0), we can evaluate the evolutionary 
effect on Sn(R;< z max ) or A^(< z m3iX ). As examples, we consider three representative CDM 
models (Table 1) whose fluctuation amplitude as is normalized so as to reproduce the abundances 
of clusters of galaxies (e.g., Kitayama & Suto 1997; Kitayama, Sasaki & Suto 1997). The results 
are displayed on a series of figures. Figure |l] shows how a is related to the comoving smoothing 
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length R in these models. Figure || displays An(R; z) as a function of z, and finally Figure || plots 
An(R;z) against R. 

The general appearance of the figures suggests that the light-cone effect is a fairly robust 
feature, although its details depend on the model. In all cases SCDM appears to give the strongest 
effect, while for OCDM and LCDM it is slightly less pronounced. Nevertheless the difference is 
fairly small and qualitatively all models behave similarly. Note also, that the magnitude of the 
correction depends on the order N, and, in accord with intuition, it is monotonically increasing for 
higher order. As expected, the light-cone effect becomes larger as z max increases, which can be seen 
in Figure |2[ Although the correction is relatively small for shallow surveys with z £ 0.2 samples, 
Atv(-R;< Zmax) becomes ^ 10% in nonlinear scales (R ~ l/i _1 Mpc). In SCDM, for instance, 
An(R',< Zmax) exceeds unity for N > 6 for the entire dynamic range plotted. Furthermore 
Figure || indicates that even if the hierarchical ansatz is correct, i.e., Sn(R',z) is independent of 
R, the light-cone effect should generate apparent scale-dependence, since the correction behaves 
differently at different scales at a given redshift. 

The future SDSS will be able to measure the moments of the galaxy density field with 
unprecedented accuracy. Unless unforeseen systematics exists, it will determine them with less 
than a few percent error for N < 3 and 10% for N = 4 between 1 and 50/i _1 Mpc (see Colombi, 
Szapudi, & Szalay 1997 for details). According to Figures ^ and [|, the light cone effect will be 
much smaller than these errors, or at most of the same order, depending on the scales and models. 
The correction could be potentially non-negligible only at the smallest scales. A similar conclusion 
can probably be drawn about the 2dF survey. On the other hands, for future deep surveys, which 
should aim at smaller scales especially if carried out by the Next-Generation Space Telescope, 
our calculations will be of utmost importance. According to Figure [3| the required correction can 
range from up to unity for S3 through factors of few for Sq to factors of hundred for Sio. 
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Research of the Ministry of Education, Science, Sports and Culture of Japan (07CE2002) to 
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Table 1: CDM model parameters 
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Fig. 1.— a(R) is plotted against log 10 i?(l/t _1 Mpc) for SCDM (solid line), OCDM (thin dotted 
line), and LCDM (thick dotted line) summarized in Table 1. 
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R = 1 h- 1 Mpc R = 10 h- 1 Mpc 




0.1 1 0.1 1 

redshift z redshift z 



Fig. 2. — log 10 An(R; z) are shown as functions of log 10 z at R = l/i _1 Mpc (/e/t panels) and 
lOh^Mpc (n^/ti panels); SCDM (top panels), OCDM (middle panels), and LCDM (bottom panels). 
The family of curves display different orders from iV = 3 ... AT = 10 monotonically upward; for 
iV = 3, and iV = 7 is plotted with thick lines for orientation. 
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Fig. 3. — logio ^n(R] z) are displayed as functions of log 10 R at z = 0.2 {left panels) and 1.0 {right 
panels); SCDM {top panels), OCDM {middle panels), and LCDM {bottom panels). The family of 
curves is the same as for Fig. 2. 



